Generating two-photon entangled states in a driven two-atom system 
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We describe a mechanism for a controlled generation of a pure Bell state with correlated atoms that involve 
two or zero excitations. The mechanism inhibits transitions into singly excited collective states of a two-atom 
system by shifting them from their unperturbed energies. The shift is accomplished by the dipole-dipole inter- 
action between the atoms. The creation of the Bell state is found to be dependent on the relaxation of the atomic 
excitation. When the relaxation is not present or can be ignored, the state of the system evolves harmonically 
between a separable to the maximally entangled state. We follow the temporal evolution of the state and find that 
the concurrence can be different from zero only in the presence of the dipole-dipole interaction. Furthermore, 
in the limit of a large dipole-dipole interaction, the concurrence reduces to that predicted for an X-state of the 
system. A general inequality is found which shows that the concurrence of an X-state system is a lower bound 
for the concurrence of the two-atom system. With the relaxation present, the general state of the system is a 
mixed state that under a strong dipole-dipole interaction reduces the system to an X-state form. We find that 
mixed states admit of lower level of entanglement, and the entanglement may occur over a finite range of time. 
A simple analytical expression is obtained for the steady-state concurrence which shows that there is a threshold 
value for the dipole-dipole interaction relative to the Rabi frequency of the driving field above which the atoms 
can be entangled over the entire time of the evolution. 

PACS numbers: 03.65.Yz, 03.67.Bg, 42.50.Dv, 42.50.Hz 



I. INTRODUCTION 

It has been known for many years that two types of max- 
imally entangled Bell states can be generated in a system 
composed of two two-level atoms; the so-called spin anti- 
correlated states that are linear superpositions of single exci- 
tation product states, and spin correlated states that are linear 
superpositions of double and zero excitation product states. 
Recent theoretical and experimental work has demonstrated 
that single excitation Bell states can be generated in a system 
of two Rydberg atoms by the dipole-dipole blockade mech- 
anism [ 1 ] . The mechanism is often referred to as Rydberg 
blockade and bears on the elimination of the simultaneous ex- 
citation of the atoms by shifting the double-excitation states of 
the system, an effect attributed to the long range dipole-dipole 
interaction characteristics of Rydberg atoms IBl- llOn . Rydberg 
atoms are highly excited atoms which have large sizes and 
therefore can have huge dipole moments, proportional to n 2 , 
where n is the principal quantum number lUlll . Because the 
dipole moments of Rydberg atoms are so large, the atoms can 
strongly interact with each other even at large distances. 

The physical origin of the dipole-dipole blockade is ex- 
plained clearly in terms of collective states of the two-atom 
system. These states provide a more natural basis for inter- 
acting atoms Il2l - [l5ll . The effects of the dipole-dipole in- 
teraction include in general the creation of maximally entan- 
gled single excitation states and the shift of the states from 
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the single-atom energy Il6ljl8ll . When the double-excitation 
states are shifted form their resonances, the two-photon exci- 
tation of the system by a resonant laser field becomes sup- 
pressed. The two-photon excitation is suppressed without 
destroying the one-photon excitation. The blockade effect 
is a simple process for creation of single-excitation entan- 
gled states, and the preparation of the entangled states via 
the dipole blockade has recently been demonstrated experi- 
mentally lfl9U22ll . The presence of only single excitations of 
the system is manifested by the photon antibunching effect, 
which signifies that at g iven time only a single photon is emit- 
ted by the system Il23l - l26ll . In the experiment of Gaetan et 
al. IU9I1 . a modification of the Rabi frequency of the driving 
laser field by \[2 has been observed that is recognised as a 
signature of not only the dipole-dipole blockade effect, but 
also of the creation of a single excitation entangled state. In 
an earlier dipole blockade experiment, Heidemann et al. ll20ll 
have observed the characteristic ^/N scaling of the Rabi os- 
cillations between the ground state and the single excitation 
multi-atom entangled state of a laser driven mesoscopic en- 
semble of TV ultracold Rydberg atoms. In the experiments of 
Zhang et al. i2lll and Wilk et al. l22tl . an entanglement be- 
tween two Rydberg atoms has been demonstrated by measur- 
ing the state fidelity of F = 0.71 and F = 0.75, respectively, 
that are well above the threshold of F = 0.5 required for 
quantum entanglement. The experiments clearly demonstrate 
that the dipole-dipole blockade can deterministically generate 
entangled states in an atomic system. 

The notion of blockade of multi-photon excitations is not 
restricted to the dipole-dipole blockade, but has been extended 
to photon blockade for the transport of light through an opti- 
cal system |[27l - l29ll . and Coulomb blockade of resonant trans- 
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port of electrons through small metallic or semiconductor de- 
vices OOll . The photon blockade mechanism, similar to the 
dipole-dipole blockade, prevents absorption of more than one 
photon by an optical system. In analogy, the Coulomb block- 
ade prevents transport of more than one electron through a 
metallic or semiconductor device. 

The dipole-dipole blockade mechanism described above 
applies to the process of the creation of maximally entan- 
gled single excitation Bell states only. It is the purpose of 
this paper, therefore, to propose a mechanism which might 
be useful in achieving entangled two-photon Bell states in a 
system composed of two coherently driven two-level atoms. 
As we have already mentioned, these states are linear super- 
position of product states with double and zero excitations. 
We shall demonstrate that the states could also be gener- 
ated with the help of the dipole-dipole interaction between 
the atoms. Specifically, we analyze the dipole-dipole inter- 
action as a blockade effect for a single photon absorption of 
the laser field. We work in the collective basis of the system, 
and find that in the limit of a large dipole-dipole interaction, 
the collective four level system reduces to an effective two- 
level two-photon system. The state of this reduced system is 
then obtained analytically and the nature of the state is fully 
analyzed. We show that the creation of a pure two-photon Bell 
state is dependent on the relaxation of the atomic excitation. 
When the relaxation is not present or can be ignored, the state 
of the system evolves harmonically between a separable to the 
maximally entangled two-photon Bell states. With the relax- 
ation present, the general state of the system is a mixed state 
that under a strong dipole-dipole interaction reduces the sys- 
tem to an X-state form. We find that mixed states admit of 
lower level of entanglement, and show that there is a thresh- 
old value for the dipole-dipole interaction relative to the Rabi 
frequency of the driving field above which the atoms can be 
entangled over the entire time of the evolution. 

The paper is organized as follows. In Sec. [TTJ we give a 
qualitative explanation of the mechanism for the creation of 
two-photon Bell states in a driven two-atom system with the 
help of the dipole-dipole interaction. A detailed calculation 
of the evolution of the system isolated from a dissipative en- 
vironment is studied in Sec. [HI] Analytical expressions are 
obtained for the probability amplitudes of the collective states 
of the system along with the expression for the concurrence. 
We show how the shift of the single excitation states may lead 
to a maximally entangled state involving only the ground and 
double excited states of the system. In Sec. IIII Bl we present 
the conditions for the system to reduce to a two-level two- 
photon system described by the density matrix in an X-state 
form. Section[lV]deals with a practically more realistic model 
in which the atoms interact with a vacuum reservoir. The in- 
teraction results in the dissipation of the atomic excitation and 
coherence, and we examine the effect of the dissipation on 
the transient evolution and stationary properties of the con- 
currence. The remarkably simple analytical expression is ob- 
tained for the steady-state concurrence. We summarize our 
results in the concluding Sec. [VI] Finally, in the Appendix, 
we present the full set of equations of motion for the density 
matrix elements and their steady-state solutions. 



II. THE MODEL 

We begin with a qualitative explanation of the concept of 
the generation of two-photon Bell states in a driven two-atom 
system with the help of the dipole-dipole interaction. Let us 
consider a system composed of two identical non-overlapping 
atoms, separated by a distance r\2 = |r*2 — f\ | and coupled to 
each other through the dipole-dipole interaction. Each atom 
is modelled as a two-level system (qubit) with the ground 
state \gi) and the excited state |e;), (i = 1, 2) separated by a 
transition frequency luq. In the absence of any external fields, 
e.g. a driving laser field, the Hamiltonian of the system is of 
the form 



H = huj (S? + S z 2 ) + H 
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where S+ = \e t )(g t \, S t = \g l )(e l \ and S- = (|ei)(e;| - 
\gi)(gi\)/2 are, respectively, the raising, lowering and energy 
difference operators of the ith atom, and fly is the dipole- 
dipole potential between the atoms. The potential depends 
on the distance between the atoms, mutual orientation of the 
atomic transition dipole moments, and the orientation of the 
dipole moments in respect to the interatomic axis. For the case 
of mutually parallel dipole moments, the potential is given 
by CM! 
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where 7 is the spontaneous emission rate of the atomic excita- 
tions, assumed the same for both atoms, 9 is the angle between 
the dipole moments of the atoms and the direction of the inter- 
atomic axis, k = u>o/c, and is the distance between the 
atoms. 

We can write the Hamiltonian dT} in a matrix form using 
basis states of two non-interacting atoms that are four product 
states 



|l) = |fli>8>|sa}, 
|3) = Isi) ® |e 2 ), 



|2) = |ei)®|e 2 ) 
|4) = |ei) ® |<?2> 



(3) 



and find 



H = h 



( -luq \ 

uj 

fii2 

V flu / 



(4) 



Notice that in the absence of the dipole-dipole interaction, the 
single excitation states |3) and |4) are degenerate in energy 
with E3 = E4 = 0, and are separated from the zero |1) and 
double 1 2) excitation states by ?wjq, as illustrated in Fig. 02a). 

In the presence of the dipole-dipole interaction, the Hamil- 
tonian (O is not diagonal, and a diagonalization results in the 
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so-called collective or Dicke states lfl2l - [l5ll 

|ei) <g> |e 2 > 
|ei> O |.92> + 8> |ca», 



|fl> = Iffi) ® |fla>, |e) 

lei) ® |<? 2 ) - |<?i) (8) |ea», 



1 ' V2 X 
with the corresponding energies 

E„ = —huiQ, E e — hujQ, E s = Ml\2, E a 



(5) 



(6) 



We see that the dipole-dipole interaction lifts the degeneracy 
of the single excitation product states and leads to two non- 
degenerate maximally entangled states \s) and \a) separated 
in energy by 2M7i2, as illustrated in Fig.[TJb). 
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FIG. 1: Schematic illustration of the idea of the creation of two- 
photon entangled states in a driven two-atom system, (a) In the ab- 
sence of the dipole-dipole interaction between the atoms, f2i2 = 0, 
the driving laser field of frequency ujl = u)o is on resonance with 
both, one and two-photon transitions of the two-atom system, (b) In 
the presence of the dipole-dipole interaction, 7^ 0, and then the 
intermediate states are shifted from the one-photon resonance leav- 
ing only the two-photon transition on resonance with the driving field 
frequency. 

We are now in a position to explain qualitatively the idea of 
the creation of two-photon entangled states with the help of 
the dipole-dipole interaction. The idea is illustrated in Fig. [TJ 
in which the schematic energy-level diagram of a two-atom 
system is shown. Independent of the basis used, a system 
composed of two two-level atoms is equivalent to a single 
four-level system with a ground state, two intermediate states 
separated from the ground state by energy Hojq and an upper 
state separated from the ground state by energy 2hojQ. As- 
sume that the two-atom system is driven by an external laser 
field of frequency tu l resonant with the atomic transition fre- 
quency, i.e. lul = ujq. Consider separately two cases, the 
absence and the presence of the dipole-dipole interaction be- 
tween the atoms. In the absence of the dipole-dipole inter- 
action, the laser field drives on resonance both, the one- and 
two-photon transitions of the system, as seen from Fig. [TJa). 
In this case, all states of the system are populated by the laser 



field. Since the one-photon coupling dominates over the two- 
photon coupling, the intermediate states |3) and |4) are more 
populated than the upper state 1 2) . The situation changes in 
the presence of the dipole-dipole interaction. The interaction 
shifts the intermediate states from the one-photon resonance, 
as shown in Fig.[TJb). As a consequence, the one-photon tran- 
sitions become off resonance to the driving field frequency, 
but the two-photon transition remains on resonance. Thus, in 
the presence of the dipole-dipole interaction, a resonant laser 
field effectively couples only to the two-photon transition of 
the two-atom system. 

A more quantitative attempt is presented below, where we 
show that the shift of the intermediate states by the dipole- 
dipole interaction can lead to the population of the upper 
state |e) without populating the intermediate states \s) and 
I a), which then may result in a pure maximally entangled two- 
photon state of the system. 



III. CREATION OF TWO-PHOTON ENTANGLED STATES 

Let us first consider the case when the atoms are isolated 
from the environment. In this case, there is no relaxation of 
the atomic excitation, but there could still exist a non-zero 
dipole-dipole coupling between the atoms. In other words, 
there are no losses due to spontaneous emission and therefore 
the evolution of the system can be determined by the evolu- 
tion of a pure state of the system. In practical terms, it could 
be done by placing the atoms inside separate cavities and to 
arrange the coupling between the cavities through an opti- 
cal waveguide Oil I32I1 . Alternatively, one could trap single 
ions or atoms inside separate potential wells and use an ad- 
ditional tr app ed ion as antennae to enhance the dipole-dipole 
coupling 113311 . 

Assume that the two-atom system described by the Hamil- 
tonian (Q~|i is subjected to a continuous driving by an exter- 
nal coherent laser field. The strength of the driving is char- 
acterized by the Rabi frequency Hq. The driving laser field 
has a traveling-wave character and is propagating in the di- 
rection orthogonal to the interatomic axis. The orthogonality 
of the driving field and the interatomic axis ensues that both 
atoms experience the same laser field amplitude and phase. 
The Hamiltonian for this case is given by 



H = hio Q (Sf + SI) + hJ2 VijSfSj- 

- ^f2 + S+)c-^^ + (Si + Sz)e iULt ] , (7) 

where ujl is the laser field frequency. The laser-atoms cou- 
pling part of the Hamiltonian retains only the terms which 
play a dominant role in the rotating-wave approximation 
(RWA). Anti-resonant terms which would make much smaller 
contributions have been omitted. 

When the relaxation of the atomic excitation is not present 
or can be neglected, the time evolution of the system is gov- 
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erned by the Schrodinger equation 

.. fc d|* a (t)> 



dt 



(8) 



where |\l/ s (i)) is the wave function of the driven two-atom 
system. 

Since the Hamiltonian (Q) depends explicitly on time, we 
make the unitary transformation 

|*(t))=exp[i(flo/ft)t]|*a(*)>, (9) 

where 

H = hw L {Sl+S z 2 ), (10) 
and obtain the Schrodinger equation for the transformed state 

_d\nt)) 



ih- 



dt 



ff|*(t)), 



(11) 



with the transformed Hamiltonian H of the form 

H = tiA(Sf + 5|) + m 12 (sf Sz + S^Si) 



hn (s+ + s+ 



sr + s 2 



(12) 



Here, A = ujq — ojl is the detuning of the laser frequency 
from the atomic resonance. 

The transformed Hamiltonian (fT2l does not depend on time, 
so the Schrodinger equation ([TO has the formal solution 



*(£)) = cxp -i(H/h)t |^(0)) 



(13) 



where 1^(0)) is the initial state of the system at t — 0. 

In order to find the time evolution of the wave function for 
a given initial state of the system, we write the Hamiltonian 
H in the basis of the collective states (f5]), and find 



H = h 



I -A o — n o \ 

A — f2 

-n -n f2i2 o 
V o o o -n 12 / 



(14) 



where f2 = fio/v^- 

Since the two-photon coherence, which is of the main in- 
terest here, attains maximal values for the laser frequency on 
resonance with the two-photon transition \g) <-> |e), we put 
A = in Eq. (fT4l i and readily find the following eigenvalues 
(energies) of the Hamiltonian 



e 1 = -n(n- ioi2 j , e 2 = n (n + Ui 12 
e 3 = o, E 4 = -no. 12 , 

with the corresponding eigenvectors (energy states) 

l^i) = J^(|e> + lfl»+v^l*>. 



(15) 
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(|e) + |<?» 



\1>a) = -^(|e) - \g)), |V4 



(16) 



where 



and 



a±= 2n> wlth 



n± = n±-n 12 , 



n = 



^0 + 4^12 



(17) 



(18) 



We see from Eq. (U~5l l that the single excitation antisymmet- 
ric state I a) is the eigenstate of the Hamiltonian with eigen- 
value —HQi2> the double excitation antisymmetric state \tjj3) 
is the eigenstate with eigenvalue zero, and the remaining states 
are superpositions of the single and double excitation states 
with non-degenerate eigenvalues. 



A. Time evolution of the concurrence 

We now turn to the problem of determining the form of 
an entangled state that could be created by the dipole-dipole 
interaction shift of the single excitation states and the degree 
of the resulting entanglement. As a measure of entanglement, 
we choose the concurrence that for a pure state |^ s (t)) of a 
two-atom system is given by 



c(t) = |<*,(t)|*.(t)>l, 



where 



and Oy^ (i — 1, 2) is the Pauli operator for the zth atom. 

For the system considered here, the state ^^(t)) is of the 
form 



!*.(*)} 



®4 2) i*:(*)>. 



(19) 



(20) 



|¥ a (i))=exp [-i(H /h)t]\*(t)), 



(21) 



where \ ^(t)) is given in Eq. ( TT3b . 

Given the state of the system, it is straightforward to cal- 
culate the concurrence. We have two alternative forms for 
the state |^ s (t)), the diagonal states ( TToT ) or the collective 
states ©. Of these alternatives, we will prefer the collective 
states as it is very often found that they form very convenient 
basis states to discuss properties of the concurrence. There- 
fore, we invert the transformation (Q~6} to find 

1 



1.9) = -7= l^i) - V® 1 1^2) - lv>3>) , 

|e) = A= (y/^T\ipi) - y/al\ip 2 ) + \ip3)) , 

\s)=^0\i> 1 )+^\i> 2 ), |a) = |V>4>, (22) 

and write the state vector of the system as 

|¥(t)) = C g {t)\g) + C e (t)\e) + C s (t)\s) + C a {t)\a), (23) 

where C n (t) is the probability amplitude of the ?ith state. 

The time evolution of the state vector of the system and so 
the concurrence depend, of course, on the initial state of the 
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system. Consider as the initial state the ground state \g(0)) = 
\go), corresponding to no excitation of the system at t = 0. In 
terms of the eigenstates (fTol l. the initial state is of the form 



lfl(0)} = 



■|Vi(0)) -y/al |^(0)) -IV's(O))) 



(24) 



From Eqs. (fT3l l and ( fT31 l. we then find that the time evolution 
of the coefficients C n (t) is of the form 



C ff (t) = - (a+e in -* + a_e" 



C e (t) 

C„(t) = 0. 



(a+e in - t + a_e- in + t 



1 fto 



-zf2 + t\ 



(25) 



Notice that the amplitude of the antisymmetric state is equal 
to zero for all times. This is due to the fact that for the con- 
figuration considered here, the laser field propagating in the 
direction perpendicular to the inter-atomic axis, the laser cou- 
ples exclusively to the symmetric state leaving the antisym- 
metric state completely decoupled from the driven states. For 
the initial state \go), the antisymmetric state is not populated 
at t = 0, and therefore the state will remain unpopulated for 
all t > 0. 

Using Eqs. ([T9T l and (|20| |, we find that the solution 
yields the following expression for the concurrence 



C(t)= \2C g (t)C e (t)-C 2 s (i) 

'fto N 
n 



-i 



— iQigi 



Ol2 

2ft 



a+e 



■urn 



a_e 



(26) 



It is easily verified that the concurrence vanishes when fti2 = 
0. Thus, crucial for entanglement is the presence of the dipole- 
dipole interaction that shifts of the single excitation states 
from their resonant values. The concurrence is different from 
zero and contains exponentials oscillating in time with three 
different frequencies, fti2, 2ft_, and 2ft + . Since, in general, 
the three frequencies are not commensurate, the evolution of 
concurrence is not periodic. However, for some values of the 
ratio fto/fti2, the evolution can be periodic. This is shown in 
Fig. |2] where we plot the concurrence versus time t for several 
different values of the ratio fto/fti2. For t = the system is 
separable, C(0) = 0, due to our choice of the initial states of 
the atoms. Immediately afterwards, the concurrence begins to 
increase and depending on the ratio fto/fti2, it oscillates pe- 
riodically or non-periodically in time. For these values of the 
ratio fto /CI12 at which the oscillations are periodical, the con- 
currence reaches the maximal value of C(t) = 1 that is seen to 
occur periodically over all times. For the case of non-periodic 
oscillations, the maxima are smaller than one. Thus, the ef- 
fect of the non-periodicity is clearly do decrease the amount 
of entanglement. 




0.2 0.4 0.6 0. 



1 1.2 1.4 



FIG. 2: Time evolution of the concurrence C(t) for flo = IO7 and 
different values of the ratio u = flo/flia: u = v/3/4 (solid line), 
u — 1.5 (dashed line), u = a/15/4 (dashed-dotted line), u = 2.5 
(dotted line). 



Let us examine the concurrence a little more closely. Con- 
sider separately the cases of periodic and non-periodic oscilla- 
tions. As we have already noticed from Fig. [2] and it can also 
be seen from Eq. (|26| |. for the case of periodic oscillations, 
the concurrence becomes unity at certain times satisfying the 
condition 

Vt\2t = 7r and Qt = rm, n = l,2, ... (27) 
which can happen only for the discrete values of the ratio 



° / 2 1 in 

— — = ( r , 71 = 1,2,. 

Oi2 V 4' 



(28) 



It is easily verified that at the times satisfying the condition 
d27i i, the probability amplitudes (l25"l l reduce to 



c g (t) = -L e (- 1 )"- 1 -/4 ; Ce{t) = __L e (-D"W4 j 



C s {t)=0, 7i = l,2,. 



(29) 



This shows that at these particular times, the system is in a 
superposition of the ground \g) and the upper |e) states with 
no excitation of the symmetric state \s). 

When Eq. ( 1291 is used in Eq. d23T l. we readily find for the 
state vector of the system of the resulting maximally entangled 
two-photon Bell state 



l*(*n)> 



; (-l) n MT/4 

V2 



|e> + (-r/%)), n = l,2, 



(30) 



We now turn to a detailed analysis of the concurrence for 
the case of non-periodic oscillations. In this case, the com- 
peting effects of one- and two-photon transitions modify the 
Rabi oscillations and cause increasing distortions of the con- 
currence. However, we can reduce the destructive effects 
of the competing transitions by taking large values of the 
dipole-dipole interaction, fii2 ^> ft. It is easily verified from 
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Eq. (l26l i. that in the limit of fii2 3> ^, the amplitudes of 
the terms oscillating with frequencies H12 and 2il + are very 
small, and then the concurrence can be approximated by 



m « - 



- 1 + exp 



29.1 



We see that at times 



12 



2Q§ 



n = 1, 3, 5, . . 



(31) 



(32) 



the concurrence is close to its maximal value of C(t) = 1. 

At these times, the state of the system reduces to the maxi- 
mally entangled two-photon Bell state 



!*(*)} = - 



«r/ 4 
75" 



(|e)-%». 



(33) 



It is seen that the system evolves between the ground \g) and 
the doubly-excited state |e) and at some discrete time can be 
found in the maximally entangled two-photon Bell state. The 
maximal entanglement is due to the two-photon coherence. 

It is straightforward to show using Eq. ((25) that the one- 
and two-photon coherences evolve in time as 



,(t)=C,(*)C;(i) = ig[e 



-iO_t _ 
2iSlt 



e iSl+t 



+a_ (c~ 2lQt - 1) - a + (e 
p es (t) = C e (t)C* s (t) = [-e-^- 4 + e^+ 4 



-2iQt 



1) - a. 



P BB (t) = c g (t)c;(t) = -U~ 



1 




sin fit 



+ —i (a_ sinf2 + t — a + sinO_t) , 



(34) 



where p ge (t) = p ge {t) exp(— 2icJot) is the slowly varying part 
of the coherence. 

It follows from Eqs. d28l i. (l32l and ( l34l that the one-photon 
coherences vanish, p gs {t) — p se (t) = 0, and the two-photon 
coherence |p ge (t)| becomes maximal, |/5 se (i)| = 1/2, at times 
the concurrence is maximal. Thus, we must conclude that the 
entanglement is created by the two-photon coherence with no 
population of the single excitation state \s). 



to zero and then the density matrix of the system takes the X- 
state form 



P(t) = 



f Pag® 

P aa(t) 








Pge(t) 









V Peg(t) 



P„(t) 

Pe*(t) 



(35) 



The concurrence of the system determined by the A-state den- 
sity matrix is readily found to be 



C x (t) = 2\p ge (t)\ - p ss (t) 

= 2\C g {t)C:{t)\-\C s (i)\\ 



(36) 



where we have taken into account that p aa {t) = 0. 

We wish to compare the concurrence C x (t) with the concur- 
rence C(t), Eq. (|26| |. and see whether the concurrence of the 
system discussed above is larger or smaller than that predicted 
by the X-state system. From Eq. d26l i. we have 



C{t)=2\C g {t)C e (t)-C 2 s {t)l 
and then applying the inequality 

\Z\Z2 - zfl > \ziz%\ - |z 3 | 2 , 



(37) 



(38) 



which holds for arbitrary complex numbers Z\, Z2 and Z3, we 
readily find that C(t) > C x (t). This inequality always holds 
true, so we may conclude that values of C x (t), which deter- 
mines concurrence for the X-state system, are lower bounds 
for the concurrence of the present system. In addition, if 
the phases of the complex amplitudes are chosen such that 
<t> g + 4>e — 20 s = 0, the inequality becomes equality. More- 
over, for the case of C s (t) = 0, that occurs at times the con- 
currence C(t) reaches the maximum C(t) = 1 value, we have 
C(t) = C x (t). This means that under the ideal two-photon 
excitation, the system behaves as an X-state system. 

The above considerations are illustrated in Fig. |3] which 
shows the concurrences C(t) and C x (t) as a function of time 
for two different values of SI12. For a small O12, the time evo- 
lution of the concurrences exhibits a modulation of different 
amplitudes, and the effect of an increasing SI12 is evident in 
the decrease of the difference between C(t) and C x (t). 



IV. MASTER EQUATION 



B. A-state as a lower bound for entanglement 

The results for entanglement created by the two-photon 
coherence bear interesting relation to an A-state entangle- 
ment rf34tl . Therefore, it is interesting to contrast the con- 
currence considered under the two-photon coherence with the 
concurrence of the A-state system. Since the dipole-dipole 
interaction shifts the single excitation states from resonance 
with the driving field, the one-photon coherences are all equal 



The above discussion of the creation of the two-photon Bell 
state describes a simplified case, which ignores dissipative ef- 
fects of spontaneous emission resulting from the coupling of 
the atoms to an external multi-mode reservoir. We now ex- 
tend the model to include the dissipation in the system. In this 
case, a general state of the system is a mixed state and the dy- 
namics of the atoms are then determined by the evolution of 
the density operator of the system. The density operator sat- 
isfies the master equation which in the interaction picture can 
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interatomic axis 



FIG. 3: Time development of the concurrences C(t) (solid line) and 
C x (t) (dashed line) for flo = IO7 and different values of fii2: (a) 
^12 = 7, and (b) Q12 = IOO7. In the case (b), both concurrences 
overlap. 



be written as 



= - iA [SI + Sl,p] - i «ii [StSj, P ] 



1 2 

-ino52[st + sr, P ] 



(=1 



E ^{pStSj + S+Sjp - 2S- P S+). 



(39) 



The term of the master equation, involving 7^, represents the 
evolution of p due to dissipation in the atomic system. The 
term contains contributions, proportional to 7^ = 7 repre- 
senting the damping of the ith atom by spontaneous emission, 
assumed to be independent of i, so that both atoms are equally 
damped by the field. Apart from the contribution of the in- 
dividual atoms, the term contains contributions, proportional 
to 7ij (i j), that represent the collective damping result- 
ing from the mutual exchange of spontaneously emitted pho- 
tons through the common reservoir lll2l - ll5ll . The parameter 
7*j (* 7^ i) depends on the distance between the atoms and 
the orientation of the atomic dipole moments in respect to the 



1 -3 cos 2 



kTij 

cos(krij) sin(fcrjj) 



(kr i:j ) (km)" 



(40) 



In general, 7^ involves terms oscillating at frequency fcry, 
multiplied by inverse powers of fcry ranging from (fcr.y ) -1 
to (krij)~ 3 . Note that 7^ is symmetric under the exchange 
of the atoms, i.e., 721 = 712- For small atomic separations, 
krij <C 1, the collective damping 7^ approaches 7, whereas 
for large separations 7^ vanishes. The later corresponds to 
the case where the atoms are independently damped by the 
reservoir. This case is equivalent to the situation where the 
atoms are damped by their own independent reservoirs. 

We now employ the master equation (l39l to find equations 
of motion for the matrix elements of the atomic density oper- 
ator p. We use the collective states basis ([5} as the represen- 
tation of the density operator. As in the previous sections, we 
assume that both atoms experience the same amplitude and 
phase of the driving field and we focus on the two-photon res- 
onance, i.e., we put A = 0. In so doing, we find that the 
equations of motion split into two independent sets, one com- 
posed of nine inhomogeneous and the other composed of six 
homogeneous coupled equations. The sets together with their 
steady-state solutions are listed in Appendix A. 



A. Effect of spontaneous emission on entanglement 

We wish to examine the time development of the concur- 
rence in the presence of spontaneous decay of the atomic exci- 
tation and coherence. In this case a general state of the system 
is a mixed state described by the density operator p. 

The concurrence of a mixed state of a two-atom system is 
defined as 0511 



C = max ( 0, v Ai — v A2 — v ^3 ~ V ^ 



(41) 



where Ai are the the eigenvalues, putted in decreasing order, 
of the matrix 



with 



P 



R = pp, 



UyP Uy ® Cy 



(42) 



(43) 



and p* denotes the complex conjugate of p. 

The concurrence is specified by the density matrix of a 
given system and thus can be determined from the knowl- 
edge of the density matrix elements. Assume that initially, 
prior the laser field was turned on at t = 0, the atoms were 
in their ground states, i.e., p gg (0) = 1 and the other matrix 
elements equal to zero. It is easily verified from Eq. (IA5t 
that in this case, the coherences between the triplet states and 
the antisymmetric state are equal to zero for all times. As a 
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consequence, the density matrix of the system, written in the 
collective basis (\g), |e), \s), \a)), takes the simplified form 



P(t) = 



/ Pgg(t) Pge(t) Pgs{t) 

Peg(t) Pee(t) Pes(t) 

Psg{t) p S e(t) Pss(t) 

V Paa {t) 



(44) 



The time evolution of the density matrix elements is found 
by solving the set of nine coupled differential equations 1A41 . 
These equations are cumbersome for an analytical solution 
because of the coupling between the populations and coher- 
ences. Therefore, we use a numerical method to find the 
evolution of the density matrix elements from which we then 
compute the time evolution of the concurrence. 

Note that the density matrix (l44t is not diagonal in the basis 
of the collective states. This means that, in general, the collec- 
tive states are not the eigenstates of the system. In principle, 
it is possible to find the diagonal states simply by direct diag- 
onalization of the matrix (l44l . However, diagonal states such 
obtained are complicated in form and thus difficult to inter- 
pret. An alternative way is to compare the general state (Pffl 
with approximate states in which the system could be found 
under the pure two-photon excitation. 

Of these alternatives, we consider a state described by an 
approximate density operator of the form 



pit) = (1 - 4p ao (t))|*(f)>(*(t)| + P aa(t% 



(45) 



where 1^(4)) is the pure state of the system, calculated in 
Sec. [nil an d I is the 4x4 identity matrix. The approxi- 
mate density operator represents a state of the system that is 
an incoherent superposition of the pure state ^(t)) and the 
antisymmetric state \a). This choice is suggested by an ob- 
servation from Eq. (IA41 > that the spontaneous decay couples 
the antisymmetric state to the triplet states. This means that 
this coupling may play an important role in the creation of the 
two-photon Bell states. 

We also compare the general state (l44t with a state de- 
scribed by the density matrix of the AT-state form. Therefore, 
we consider the criterion C x for entanglement of an A-state 
system given by ll36Tl 

C x {t) = max{0,2\p ge (t)\ 

-y/\pss(t) + p aa (t)] 2 - [2Rep sa (i)] 2 } . (46) 

Since for the general state p sa (t) = 0, the criterion (|46*] | sim- 
plifies to 

C x (t)=m&x{0,2\p ge (t)\-[p ss (t)+p aa (t)}}. (47) 

We now present some numerical calculations that illustrate 
the effect of spontaneous emission on the pure two-photon 
Bell state created with the help of the dipole-dipole interac- 
tion. The time development of the concurrence is calculated 
for the actual state described by the density matrix (PPfl and, 
for comparison, for an approximate state described by the den- 
sity matrix (PES! and also for an X state of the system. 
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FIG. 4: Time development of the concurrence for three different 
states of the system. The solid line shows the concurrence of the 
pure state of the system decoupled from the reservoir. The dashed 
line is the concurrence of the actual state of the atoms coupled to 
a common reservoir, and the dashed-doted line is the concurrence of 
the approximate st ate of the system, Eq. d45l l. In each case flo = IO7 
and f2i2 = fio In last two cases, the distance between the 
atoms is ria = 0.078A and the dipole moments are polarized in the 
direction perpendicular to the interatomic axis, 6 = tt/2. This choice 
of the parameters has been made to get = 11. 487 = flo/\/3/4. 



Figure|4]shows the time development of the concurrence of 
the system of two atoms located at a small distance, ri2 = 
0.078A, and coupled to a common reservoir. Also shown are 
concurrences of the system decoupled from the reservoir and 
that determined by the density matrix (l45t . It is easy seen that 
the concurrences evolve in decidedly different ways. In partic- 
ular, the concurrence of the state approximated by the density 
matrix d45l l does not reproduce the concurrence of the actual 
state of the system. It is rather close to the concurrence of the 
pure state of the system of atoms decoupled from the reservoir. 
The concurrence oscillates in time with the same frequency as 
that of the isolated system and persists over many Rabi pe- 
riods. The amplitude of the oscillations decreases slowly in 
time and is damped out in a time of order (7 — 7i2) _1 - 

The reason for this feature of the approximate state of the 
system can be understood by referring to the asymmetry in- 
troduced to the system by the competing effects of sponta- 
neous emission from the upper state |e) to the symmetric and 
antisymmetric states. It is easily verified from Eq. (lA4t that 
the antisymmetric state \a) is populated by spontaneous emis- 
sion from the state |e) at a sub-radiant rate 7 — 712, which for 
small distances between the atoms is much smaller that the 
superradiant rate, 7 + 712, at which the symmetric state \s) 
is populated. This means that the redistribution of the noise 
in the system is not isotropic. Consequently, the symmetric 
and antisymmetric states are affected differently by sponta- 
neous emission. A considerable part of the vacuum noise is 
transferred to the symmetric state with only a small fraction 
being transferred to the antisymmetric state. Thus, the popu- 
lation p aa (t) is not a good measure of the noise distribution 
between the states of the system of two atoms coupled to a 
common reservoir. 
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FIG. 5: The time evolution of the concurrence of the actual state 
of the system C(t) (solid line) and the approximate X state C x (t) 
(dashed line) for the same parameters as in Fig. |4] 

The actual mixed state of the system is much better approx- 
imated by the X state density matrix. It is convincingly seen 
from Fig. [5] where we graph the concurrence C(t) of the actual 
state and the concurrence C x (t) of the X state of the system 
for the same parameters as in Fig. |4] The concurrences behave 
in a similar fashion that C x (t) remains quite close to C(t) for 
all times. The entanglement lives over a restricted time range 
with no oscillations present even though the Rabi frequency is 
high. This is an example of the phenomenon of sudden death 
of entanglement 07H39ll . Notice that C x (t) does not exceed 
the actual concurrence C(t), same feature as predicted above 
for the pure state, see Sec. IIIIBI 



quite well the actual state in the case of independent reser- 
voirs, i.e., when 712 = 0. We illustrate this feature in Fig. [6] 
where we plot the time evolution of the concurrence for two 
different arrangements of the coupling of the atoms to an ex- 
ternal reservoir. It is seen from Fig. |6]that the concurrence 
of the approximate state d45l l coincides quite well with the 
concurrence of the actual state of the system. The reason for 
this similarity is that, in the limit of 712 = 0, the symmetric 
and antisymmetric states are equally populated by the spon- 
taneous emission from the upper state |e). This results in an 
equal redistribution of the noise between the symmetric and 
antisymmetric states. In other words, the antisymmetric state 
is now fully participating in the dynamics of the system, and 
the actual state of the system of atoms coupled to indepen- 
dent reservoirs can be well approximated by the density ma- 
trix (1451 

We close this part of the section by pointing out that fol- 
lowing the results presented in Figs. |4]|6]one would conclude 
that in the presence of spontaneous emission, the entangle- 
ment created by the two-photon coherence may exist only over 
a short time of the evolution of the system. However, a quick 
look at the figures reveals that the time range over which the 
entanglement exists depends on the value of the parameters 
involved. A question then arises, which of the parameters are 
crucial for the entanglement to survive over a long time, pre- 
sumably until steady state? In order to examine this point 
more closely, we consider the steady state solution for the den- 
sity matrix elements from which we then infer conditions for 
a non-zero steady state entanglement. 




FIG. 6: The time evolution of the concurrence of the dipole-dipole 
interacting atoms for two different arrangements of the coupling to 
an external reservoir. The solid line is the concurrence of the system 
decoupled from the reservoir. The dashed line is the concurrence 
of the actual state of the atoms coupled to independent reservoirs 
(712 = 0), and the dashed-dotted line is the concurrence of the ap- 
proximate state of the atoms, Eq. d45t . and also coupled to indepen- 
dent reservoirs. In each case, Slo = IO7 and 0,12 = Slo / \/3/4. 



Although the density matrix (1451 1 differs considerably from 
the density matrix of the actual state of the system when the 
atoms interact with a common reservoir, it does approximate 



B. Steady-state entanglement 

We now proceed to use the steady state solutions (IA71 > in 
order to evaluate the concurrence of the stationary state of the 
system. To evaluate the concurrence, defined in Eq. (1411 . we 
must first find the eigenvalues of the matrix R. By taking the 
explicit form of the matrix R and making use of the steady- 
state solutions for the density matrix elements, Eq. (lA7t , we 
readily find the following expressions for the square roots of 
the required eigenvalues 



Q2 Jfi4 + 4 7 2|^ 12 |2 ± 2 7 |J7 12 | 



o 4 + 7 2 2f> 2 + ( 7 + 112 y + n 2 12 



Or 



4 04 



4 + 7 2 2ft 2 + (7+7 12 r+r! 2 



(48) 



where U12 = 712 + ^12 describes the strength of the interac- 
tion between the atoms. 

Although the above formulas look complicated, a straight- 
forward but lengthy calculation leads to a remarkably simple 
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analytical expression for the steady-state concurrence 



C(oo) = max< 0, — 



n 2 


( 7 |cr 12 | - §n») 




4 + 7 2 


20 2 + ( 7 + 7l2 ) 2 





(49) 

Equation (l49b is a general formula for the steady state con- 
currence valid for an arbitrary Rabi frequency, an arbitrary 
distance between the atoms, and for a common (712 7^ 0) or 
separate (712 = 0) reservoirs. It is seen that there is a thresh- 
old for the Rabi frequency below which the atoms are entan- 
gled in the steady state. The threshold depends only on the 
relation between the Rabi frequency and the strength of the 
interaction between the atoms. The interaction is determined 
by the collective parameters CI12 and 712. 

In the special case of independent reservoirs, i.e., when 
712 = 0, the concurrence reduces to the result recently ob- 
tained by Li and Paraoanu i40ll 



C(oo) 



0, 



(7^1 



i« 2 



(0 2 + 7 2 ) 2 +7 2 2 2 



(50) 



which also exhibits a threshold for the steady-state entangle- 
ment. The existence of a threshold was also noted by Macovei 
et al. [413, who studied numerically the stationary pairwise 
entanglement in a system composed of N atoms confined to a 
region much smaller than the resonant wavelength, i.e., in the 
limit of 712 = 7 and p s aa = 0. 

The steady-state concurrence ( |49l depends on the value of 
the Rabi frequency of the driving field relative to strength of 
the interatomic interactions. Obviously, the stationary entan- 
glement is zero for independent atoms, i.e., 712 = and 
SI12 = 0. However, as soon as 712 7^ and/or 0x2 7^ 0, 
the atoms can be entangled in the steady state. The sufficient 
condition for a steady state entanglement is to maintain the 
strength of the interatomic interactions 7 1 J7i2 1 > f2 2 /2. Thus, 
by changing the value of Cl, we may dynamically switch on or 
switch off the steady state entanglement. 

Looking at the steady-state solutions for the density ma- 
trix elements, Eq. (IA7t . we see that in the limit of a strong 
dipole-dipole interaction, Q12 ^ 7, the two-photon coher- 
ence p ge approaches a large non-zero value, p ge (oo) « — i/5 
with all other coherences being vanishingly small. Evidently, 
the source of the two-photon coherence is in the dipole-dipole 
interaction that shifts the single excitation states from the res- 
onance with the driving laser field. The strong correlations 
are reflected in the stationary state of the atoms which de- 
spite their interaction with a dissipative reservoir, decay to a 
strongly correlated state. 

The existence of the threshold for the steady-state entan- 
glement provides a clear explanation of why for the exam- 
ples illustrated in Figs. |4]|6] the entanglement was confined 
to short times. Simply, the chosen values of the parameters 
corresponded to the below threshold situations of the inter- 
atomic interaction strength 7 |i7i2| < 2 /2. If instead, we 
choose the values of the parameters such that 7IC/12I > " /2, 




FIG. 7: The time evolution of the concurrences C(t) (solid line) and 
C x (t) (dashed line) for fi = IO7, = tt/2 and r 12 = 0.06A, 
(712 = 0.97 7 , fiiz = 26.22 7 ). 



i.e. above the threshold, we then should observe an entan- 
glement for all times t > 0. This is illustrated in Fig. [7] 
where we plot the concurrences C(t) and C x (t) for the case 
of 7IC/12I > 2 /2. The concurrence builds up in an oscil- 
latory fashion and then approaches a non-zero steady-state 
value C(oo) « 0.2. Clearly, the entanglement is present for all 
times t > 0. It is also seen that the concurrence of the X state 
admits of lower level of entanglement than the actual state of 
the system, i.e., C x (i) < C(t) for all times. 



V. POTENTIAL EXPERIMENTAL SCHEME 

Finally, we discuss experimental prospects for the realisa- 
tion of the two-photon Bell states in a system of two driven 
atoms. A potential experimental system might be based on 
the experimental setup of Gaetan et al. H, O and the same 
set of parameters employed in the observation of the Rydberg 
blockade effect. The only difference would be in the tuning 
of the excitation laser field to the energy levels of the system. 
Namely, instead of satisfying the condition of zeroing the de- 
tuning between the ground state \g) and the single-excitation 
entangled state \ s), the frequency of the excitation laser should 
be kept on the two-photon resonance between the ground 
state I g) and one of the upper states shifted by the dipole- 
dipole interaction potential AE± = ±C3/rf 2 , where Ca is a 
constant depending on the geometry of the experiment |pl [l9ll . 
In this case, the frequency shift Alj± = AE±/h appears as 
a detuning of the laser field from the one-photon resonance 
\g) — > \s). In other words, a frequency shift Aw± > 7 would 
prevent the excitation of the state \s). A practical value of 
the dipole-dipole frequency shift Alj± « 50 MHz and a typi- 
cal spontaneous emission rate of highly excited states of 87 Rb 
atoms of 7/27T ~ 6 MHz satisfy the requirement of Au± > 7. 
One could argue that such a scheme differs from that consid- 
ered in this paper and illustrated in Fig. Q] However, these two 
schemes are mathematically completely equivalent to each 
other. It is seen from Eq. (IA4t that the dipole-dipole inter- 
action appears in the equations of motions for the coherences 
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p gs and p se as a detuning of the laser field from the resonance 
frequency of the transition \g) —> \s). The detuning SI12 can 
be achieved either through a shift of the entangled state |s) 
from the atomic resonance by H12 or through a shifting of the 
laser field frequency from the atomic resonance by the amount 

of SI12. 

In closing, we briefly comment about a possible experimen- 
tal observation of the signature of entanglement created by 
the two-photon coherence. We have seen that the entangled 
state created by the shift of the single excitation states involves 
states with zero and double excitations. These states are intrin- 
sically connected to correlated states known in the literature as 
pairwise atomic squeezed states ll42l - l44ll . Similar states occur 
at long times for a two-atom system in a squeezed reservoir. 
Therefore, a signature of such a state should be seen in squeez- 
ing of the emitted fluorescence field. Hence, the presence of 
the entangled states could be detected simply by observing 
squeezed fluctuations in one of the quadratures of the emitted 
fluorescence field. The fluctuations are directly measurable in 
schemes involving homodyne or heterodyne detection. 




FIG. 8: The time variation of the fidelity F(t) for fi = IO7, 9 = n/2 
and 7-12 = 0.06A, (712 = 0.977, ^12 = 26.22 7 ). 

Alternatively, one could measure fidelity of the entangled 
state created by the two-photon coherence in the presence of 
spontaneous emission. It could be done in experiments similar 
to that of Refs. I2lll22tl . where the fidelity of entangled single 
excitation states was observed. Fidelity of the actual mixed 
state of the system, determined by the density matrix (PPfl i, is 
given by 

F(t) = (*|p(t)|*> = ±[p ee (t)+p ag (t)}+Mpe 9 (t)}, (5D 

where |\&) is the pure two-photon maximally entangled Bell 
state, Eq. d33i l, created by the two-photon coherence in the 
absence of spontaneous emission. 

Figure |8] shows a plot of the fidelity versus the normalised 
time jt given by this equation for the same parameters as in 
Fig- El We see that for all times, the fidelity of the actual state 
is smaller than one, but the essential point is that the concur- 
rence exceeds the value of 0.5 required for quantum entangle- 
ment. The fidelity is greatest in the transient regime, where it 



reaches the maximum value F(t) w 0.8, and then decreases 
to its steady-state value F(t) « 0.6. Therefore, the entangle- 
ment created by the two-photon coherence is not washed out 
by spontaneous emission and should be observable in practice. 



VI. CONCLUSIONS 

We have presented a mechanism for a controlled generation 
of pure or mixed Bell states with correlated atoms that involve 
double or zero excitations. The mechanism inhibits excitation 
of singly excited collective states of a two-atom system by 
shifting the states from the one-photon resonance. In particu- 
lar, we have shown that the shift of the energy levels can lead 
the system to evolve into a pure entangled Bell state that in- 
volves two-atom states with double or zero excitations. The 
crucial for the occurrence of the entangled state is the presence 
of a non-zero two-photon coherence. The degree of entangle- 
ment and the purity of the state depend on the relaxation of the 
atomic excitation. In the absence of the atomic relaxation, the 
state of the system evolves harmonically between a separable 
to the maximally entangled Bell state. We have found that the 
concurrence can be different from zero only in the presence 
of the dipole-dipole interaction. By going into the limit of a 
large dipole-dipole interaction, we have shown that the con- 
currence reduces to that predicted for an X-state of the sys- 
tem. Furthermore, we have demonstrated that the concurrence 
of an AT-state system is a lower bound for the concurrence of 
the two-atom system. In the presence of the relaxation, the 
general state of the system is a mixed state that under a strong 
dipole-dipole interaction reduces to an A-state form. We have 
found that mixed states admit of lower level of entanglement, 
and the entanglement may occur over a finite range of time. 
The time range for the entanglement depends on the relation 
between the dipole-dipole interaction and the Rabi frequency 
of the laser field. We have calculated the steady state concur- 
rence and have found there is a threshold value for the dipole- 
dipole interaction relative to the Rabi frequency above which 
the atoms can be entangled for all times. 
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Appendix A 

The matrix elements of the atomic density operator p writ- 
ten in the basis of the collective states of the system satisfy 
two independent sets of differential equations. This is a con- 
sequence of using the collective basis, and of the fact of as- 
suming that both atoms experience the same amplitude and 
phase of the driving field. The two sets can be written in com- 
pact matrix forms as 



X = MX + 1, Y = QY, 



(Al) 
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where the vectors X and Y are of the form 



X (Pee, Pss, Paa, Pge, Peg: Pes, Pse, Pgs, Psg) ) 

Y — (Paei Pea, Pga, Pag, Psa, Pas) , (A2) 



the vector / of the inhomogeneous terms is of the form 

T= (o,o,o,o,o,o,o,-^,^) T , 



(A3) 



and M and Q are, respectively, 9x9 and 6x6 matrices of 
the coefficients of the differential equations. 

The equations of motion for the nine components of the 
vector X are 

Pee = - ^IPee + i^{p S e - Pes), 
p ss = - 2(7 + 7l2)(/0.ss - Pee) 

+ i£l(pes - Pse + Pgs ~ Psg), 



Paa — 


-2(j-J 12 ){paa- Pee), 




Pge = 


~ 27 Pge + iQ(pse - Pgs), 




Peg = 


~ 27 Peg - iQ(Pes ~ Psg), 




Pes = 


~ (37 + 712 - i£ll 2 )Pes + i&(Pss ~ 


Pee Peg), 


Pse = 


~ (37 + 712 + i^l\ 2 )p S e - i&(Pss ~ 


Pee ~ Pge), 


Pgs = 


- iVt - (7 + 712 ~ i£ll2)Pgs + 2(7 


f Jl2) Pse 



+ ifl(2p ss + p aa + p ee - p ge ), 

iD. - (7 + 712 + iftl2)Psg + 2(7 + Jl2)pe 

- iVl(2p ss + p aa + p ee - p eg ), 



(A4) 



and the equations of motion for the six components of the 
vector Y are 

Pae = -(37 - 712 - i£l\2)Pae - i^Pas, 
Pea = -(37 - 712 + i£ll2)Pea + i&Psa, 

Pga = -(7 - 712 + iSll2)Pga - 2(7 - ll 2 )p a e + i&Psa, 
Pag = -(7 - 712 - iSll2)Pag - 2(7 - Jl 2 )p e a ~ l^Pas, 
p sa = -2(7 + lVL 12 )psa + i&{Pea + Pga), 

p as = -2(7 - iVt 12 )Pas - i£l(Pae + Pag), (A5) 



where Ci = il / V%- 

The matrices M and Q are non-singular, so we can read- 
ily solve the equations (IA1I) by matrix inversion. A formal 
integration gives 



X(t) = X(0) exp(Mi) + [exp(Mi) - 1] I, 

Y(t) = Y(0) exp(Qf), (A6) 



where X(0) = X(t = 0) and F(0) = Y(t = 0) are vectors 
of initial values of the matrix elements. 

It is easy to see from Eqs. (lA4t and (IA51 > that only the com- 
ponents of the vector X can have nonzero steady-state so- 
lutions. All the components of the vector Y are zero in the 
steady-state. 

After straightforward but quite tedious calculations, we find 
that the steady-state solution for the components of the vec- 
tor X are 



Paa a r\ 



Pge 

Pis 



{Peg) 
(PLY 



P 



If) 4 

4~D 

! (o 2 + 4 7 2 ) n 2 

4 D ' 

1 7 (7 + U12) » 2 

2 D 
i 7ft 3 
2~D~' 

• tQ [n 2 + 27 ( 7 + u 12 ) 

~2~ 



D 



(A7) 



where the superscript " s" stands for the steady-state value, 



L> = rV+ 7 2 2fi 2 + ( 7 + 7i2) 2 + fi 2 2 , 



and Ui 2 = 712 + iQi 2 . 



(A8) 
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